Determining Robust Impacts of Land-Use-Induced Land Cover Changes on Surface Climate over North America and Eurasia: Results from the First Set of LUCID Experiments

The project Land-Use and Climate, Identiﬁcation of Robust Impacts (LUCID) was conceived to address the robustness of biogeophysical impacts of historical land use–land cover change (LULCC). LUCID used seven atmosphere–land models with a common experimental design to explore those impacts of LULCC that are robust and consistent across the climate models. The biogeophysical impacts of LULCC were also compared to the impact of elevated greenhouse gases and resulting changes in sea surface temperatures and sea ice extent (CO2SST). Focusing the analysis on Eurasia and North America, this study shows that for a number of variables LULCC has an impact of similar magnitude but of an opposite sign, to increased greenhouse gases and warmer oceans. However, the variability among the individual models’ response to LULCC is larger than that found from the increase in CO2SST. The results of the study show that although the dispersion amongthe models’responseto LULCC islarge,thereare a numberof robustcommonfeatures shared by all models: the amount of available energy used for turbulent ﬂuxes is consistent between the models and the changes in response to LULCC depend almost linearly on the amount of trees removed. However, less encouraging is the conclusion that there is no consistency among the various models regarding how LULCC affects the partitioning of available energy between latent and sensible heat ﬂuxes at a speciﬁc time. The results therefore highlight the urgent need to evaluate land surface models more thoroughly, particularly how they respond to a perturbation in addition to how they simulate an observed average state.


Introduction
Land use-land cover change (LULCC) via deforestation, or via conversion of natural grasslands, occurs principally for urbanization and agriculture.It is a process that probably began with human's systematic use of fire ;400 000 yr ago (Williams 2003).There is no doubt that LULCC has been geographically extensive (Defries et al. 1995;Ramankutty and Foley 1999;Klein Goldewijk 2001;Hurtt et al. 2006;Pongratz et al. 2008;Klein Goldewijk et al. 2011;Kaplan et al. 2011).LULCC affects the nature of the land surface in ways strongly determined by the type of change.The impact of removing forests for agriculture depends on the spatial scale of the change.It also depends on the type of natural vegetation removed, since a coniferous forest interacts with the atmosphere differently than a deciduous forest, and both differ significantly from a native grassland or savanna system.The impact of LULCC also depends on the type of crops implemented, whether the crops are irrigated, and whether the croplands revert back to something similar to the original vegetation over time or are maintained as crops or are transformed to urban landscapes.Large-scale LULCC affects the mean regional climate (e.g., Bonan 1997;Gallo et al. 1999;Zhou et al. 2004;Oleson et al. 2004;Hale et al. 2006;Lobell and Bonfils 2008;Mahmood et al. 2008;Fall et al. 2010) but also the likelihood of extremes (Zhao and Pitman 2002;Deo et al. 2009;Teuling et al. 2010), carbon and other trace gas emissions (Denman et al. 2007), fluxes of biologically active volatile organic compounds (Arneth et al. 2010), and the direct heating of the atmosphere through anthropogenic heat input (McCarthy et al. 2010).LULCC has also had a substantial biogeochemical effect on global climate through emission of CO 2 and other greenhouse gases (GHGs), such as CH 4 and N 2 O (Denman et al. 2007), although such impacts and exchanges are not within the scope of this paper.
The biogeochemical effect of historical LULCC on climate is reasonably well established.LULCC releases CO 2 , and therefore increases radiative forcing, although the magnitude of this increase still requires accurate quantification.In contrast, the impact of the biogeophysical effects of LULCC is more uncertain.While the land surface modeling community would basically agree that LULCC is a significant driver of climate through physical effects, there is no consensus on what this statement means.The community would agree that LULCC affects the albedo and through the albedo the energy to drive the surface energy balance.LULCC also affects how that energy is partitioned into sensible and latent heat fluxes, the turbulent energy fluxes that transfer water and heat into the atmosphere.There is a general consensus that provided LULCC is spatially coherent and on a spatial scale that is large enough, this would change the regional-scale climate significantly.However, the community could not define the scale of LULCC required to be ''large enough,'' how large the resulting change in the regional climate might be expected to be, or how the nature of the existing climate over a region might suppress or amplify the initial impacts of LULCC.Most contentious is the unresolved question of whether an LULCC-induced regional climate change teleconnects to trigger significant remote (in space) changes in climate.Some climate modeling evidence suggests that it can (Henderson-Sellers et al. 1993;Zhang et al. 1996;Gedney and Valdes 2000;Werth andAvissar 2002, 2005) and other climate modeling evidence suggests it cannot (Findell et al. 2007(Findell et al. , 2009;;Pitman et al. 2009), and this is not a problem easily resolved using observational evidence.
Another complication arises from the type of climate model used in the LULCC experiments.Experiments with prescribed sea surface temperatures (SSTs) and sea ice cover could underestimate the global effects of LULCC because air temperature over ocean, which covers 71% of the planet, is constrained to the observed state.
The global-scale impact of historical LULCC on radiative forcing through land surface albedo changes estimated as 20.2 60.2 W m 22 is small relative to the CO 2 radiative forcing of 1.66 60.17 W m 22 (Forster et al. 2007).However, this does not imply LULCC has no climatic impact, as Pielke et al. (2002) and Davin et al. (2007) have argued.A globally integrated measure of the impact of a forcing is important where a forcing has a global signature.Even if most experiments seem to show that LULCC has a negligible global signature, intense LULCC-where it has transformed large regions of the earth's surface-is spatially organized into distinct regions (Fig. 1).These include (but are not restricted to) North America, Europe, India, China, Russia, Japan, and Indonesia (all account for about 64% of all arable land on about 50% of the total land-where 60% of the global population lives1 ).So, the appropriate question is not whether LULCC has a globally averaged significant impact, rather it is whether LULCC has an impact on regions that have undergone intensive LULCC that is worth accounting for when exploring the impact of other human forcings on regional climate.
To begin to systematically address this question, the project Land-Use and Climate, Identification of Robust Impacts (LUCID) was conceived under the auspices of the International Geosphere-Biosphere Programme (IGBP)-Integrated Land Ecosystem-Atmosphere Processes Study (iLEAPS) and the Global Energy and Water Cycle Experiment-Global Land Atmosphere System Study (GEWEX-GLASS).The goal of LUCID is to determine the scale of impact of historical LULCC at the regional and global scales.A key objective is to determine those impacts that are real and robust both above the level of natural variability in each model and common across the set of seven models used in the project (see section 2).
The preliminary global-scale results from LUCID were presented by Pitman et al. (2009), who mainly focused on the Northern Hemisphere summer season.The key result was a statistically significant impact of LULCC on the simulated latent heat flux and air temperature over the regions where LULCC was changed.However, the direction of summer temperature change was inconsistent among the models, which undermines our confidence in the robustness of this result.In the case of rainfall, Seneviratne et al. (2010) noted that results from LULCC were less convincing.Four of the coupled atmosphereland models used in LUCID show a strongly significant impact on rainfall over regions of LULCC, while three models do not show impacts above the level expected by chance.Therefore, in terms of the results from LUCID, the model results do not simulate a consistent LULCC signal in terms of regional rainfall.Pitman et al. (2009) also failed to identify statistically significant impacts of LULCC remote from the actual LULCC for latent heat flux, temperature, or rainfall.Note that since LUCID did not use climate models with coupled oceans, this does not resolve the question of whether these teleconnections exist, as the coupled ocean can amplify or damp an initial regional perturbation (Davin and de Noblet-Ducoudre 2010).
A critical result from LUCID was the disagreement between the models in their response to LULCC in the simulation of summer temperature and latent heat flux change.Pitman et al. (2009) did not explore the reasons behind these differences.The goal of this paper is to provide a detailed examination of why the land surface models (LSMs) diverge in their response to LULCC and to provide advice on how to better approach realistic responses of climate models to such perturbations in the future.This paper highlights in detail where the models agree and where they disagree, and it explores the reasons for the disagreement.
Section 2 summarizes the main features of the atmosphere and land surface models used, describes the experimental protocol, the imposed vegetation changes, and the way they were implemented in individual LSMs.Section 3 compares the climatic impacts of LULCC to those resulting from the changes in sea surface temperatures, sea ice extent, and increased greenhouse gases.Section 4 discusses the resulting changes in land surface properties in each individual model, as well as their simulated impacts on fluxes and surface temperature.Section 5 focuses on the way the various models partition the available surface radiative energy into latent and sensible heat fluxes, and on how this partitioning is modified when the land cover distribution is changed.The results are discussed and summarized in section 6.

a. Models description
Seven coupled atmosphere-LSMs have been used to carry out a common set of simulations using the same experimental design.The models and their references are listed in Table 1.

b. Experimental design
All seven models undertook two sets of two simulations (Table 2), spanning a matrix of present-day (PD) and preindustrial (PI) GHG concentrations-SSTs and present-day and preindustrial land cover.In these experiments the models are forced with two different vegetation distributions (representative of 1870 or 1992).Each model carried out at least five independent simulations for each experiment to increase the capacity to determine those changes that were robust from those that reflected internal model variability.Yellow and red are used when the extent of anthropogenic areas have increased since PI times, while blue refers to abandoned lands.The two contours that are drawn on the map highlight the regions that will further be used for specific analysis (North America and Eurasia).
Unauthenticated | Downloaded 06/11/21 07:46 AM UTC From these simulations we constructed responses to LULCC by averaging the results for land cover change experiments under both present-day and preindustrial GHG concentrations/SSTs [i.e., 0.5(PD 2 PDv 1 PIv 2 PI)]; values computed in this way are termed LULCC in the following.Similarly, the average response to GHG and SST changes is evaluated as 0.5(PD 2 PIv 1 PDv 2 PI) and labeled ''CO2SST'' (the combined CO2 and SST impacts) in the following.Results from all individual members of every participating model have been averaged in all the results presented below.
Both SSTs and sea ice extent were prescribed to vary interannually and seasonally using the Climate of the Twentieth Century Project specifications [see Met Office Hadley Centre Sea Ice and Sea Surface Temperature version 1.1 (HadISST1.1)at ftp://www.iges.org/pub/kinter/c20c/HadISST/].For the vegetation distribution, each model was provided the same distribution of crop and pasture (Fig. 1), at a resolution of 0.58 3 0.58.These were constructed by Ramankutty and Foley (1999) and combined with pasture area from Klein Goldewijk (2001).The fields are similar to the ones that are now being used for the Intergovernmental Panel on Climate Change (IPCC) Coupled Model Intercomparison Project phase 5 (CMIP5) historical simulations (Hurtt et al. 2011).

c. Resulting vegetation changes
The datasets provided to each modeling group did not specify the distribution of natural vegetation (land cover types that are neither crops nor pasture).The distribution of natural vegetation had to be obtained from the datasets currently used by each LSM from their ''background'' land cover map (BLCm) or from prognostic simulations of vegetation distribution.This was unavoidable, as climate models and their LSMs are commonly developed together with an assumed natural vegetation map.It is simply not feasible to require all modeling groups to use a common natural vegetation map, since it would require a restructuring of the LSM and a full-scale reevaluation of the coupled model.
As an implication, each LSM therefore integrated the crops and pasture distributions provided for both periods into its own land cover map.This has resulted in rather different land cover distributions (Fig. 2a) for  (TESSEL).Moreover, although the areas covered with crops increase in all models from preindustrial times to the present (Fig. 2b), this increase varies from about 30% in CCSM-CLM, which is less than the 35% originally prescribed, to about 55% in CCAM-CABLE.
The models therefore all experience temperate deforestation but at varying degrees.For example, for North America, temperate deforestation ranges from 12% in ECHAM5-Jena Scheme for Biosphere-Atmosphere Coupling in Hamburg (JSBACH) to 40% in EC-EARTH-TESSEL (Table 4).In EC-EARTH-TESSEL, Simplified Parameterizations, Primitive-Equation Dynamics (SPEEDY)-Lund-Potsdam-Jena Model for managed land (LPJmL), ARPEGE-ISBA, and L'Institut  was .51% in a grid cell, then the grid cell is assigned as crop.If the grid cell is crop in BLCm but not according to the crop and pasture map, then the grid cell is changed to the nearest natural vegetation.The aim is to preserve the differences in the crop and pasture extent in the PI and PD land cover maps.
Pierre-Simon Laplace (IPSL)-Organizing Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE), crops have expanded predominantly at the expense of forests, while in CCAM-CABLE, CCSM-CLM, and ECHAM5-JSBACH, the reduction in herbaceous areas is as large or even larger than the reduction in forests.
For the LUCID experiments, five out of seven models have chosen to proportionally reduce/extend all natural vegetation types on the part of the grid cell that is not occupied by crops and/or pasture, while the remaining two have first reduced the grassland area and then forests.Other differences between the resulting maps come from 1) the number of vegetation types accounted for per grid cell and 2) the BLCm.Regarding the first point, CCAM-CABLE, for example, only considers the dominant vegetation type per grid cell.If one type exceeds 51%, then it is set to occupy the whole grid cell.This is why the increase in crop area is much larger in this model than in all others.CCAM-CABLE, moreover, groups both pasture and crops into a single vegetation type and both increase from preindustrial times to present day.The extent of herbaceous types decreases in most other models between the two periods, mainly because natural grassland and pasture are most often grouped within one vegetation type.Since natural grassland is severely reduced to allow for the growth of crops, the total grassland area is decreased.As with the choice of whether to grow crops at the expense of forests or herbaceous plants, there is no known ''right'' choice in the number of vegetation types per grid cell and therefore there is no right choice in how to implement LULCC in a climate model.However, our findings of large differences in the deforestation rate, for example, suggest that we should, collectively, focus on appropriate strategies to implement LULCC in our models.The remaining differences between the resulting maps come from the various configurations of the BLCm among the modeling groups.Various examples exist, including the following: 1) The LSM uses a potential vegetation map, that is, one that only includes natural vegetation types.The crop and pasture extent will therefore, for each individual year, replace some or all natural types (e.g., CCSM-CLM, ECHAM5-JSBACH).
2) The LSM uses an observed present-day land cover map that includes both natural and anthropogenic vegetation types.Reconstructing the historical land cover map means deciding which natural vegetation types, crops and pasture would be replaced or deciding what types existed prior to the current distribution of anthropogenic vegetation.Moreover, the distribution of crops and pasture in the BLCm may be different from the present-day distribution provided We emphasize that the way each modeling group has implemented LULCC is reasonable at this time because the land surface modeling community has not recognized before how many ways there are to implement LULCC, and the required research has not yet been done to establish a preferred methodology.However, while each approach is reasonable, they are different and those differences inevitably have an impact on how LULCC affects the near-surface climate.
3. How does the models' response to the LULCC signal compare to the response to other changes?
The four sets of simulations described in section 2b allow a comparison of the sensitivity of models to both LULCC and the combined changes in atmospheric greenhouse gases, and CO2SST.To do this, we use the framework of the surface energy balance equation that is implemented in all LSMs.In the absence of snow and on longer time scales, the basic form of this equation is where all fluxes are in watts per square meter, Q S is the shortwave radiation incident at the land surface, a is the surface albedo, Q Ld is the downwelling infrared radiation, Q Lu is upwelling infrared radiation, Q H is the sensible heat flux, Q E is the latent heat flux, and Q G is the flux of heat being transmitted to deeper soil layers.The left-hand side of this equation is also known as net radiation.In our analysis we use available energy (Q A ; and not net radiation), defined as (2) Our choice for not including Q Lu in the available energy (as is more traditionally done) is because we wish to separate Q A into turbulent fluxes on the one hand and long-wave cooling on the other.We also define the total amount of energy exchanged as a turbulent energy flux (Q T ) as The responses of Q A and ambient air temperature to both sets of changes are displayed in Fig. 3 for North America and Eurasia for all seasons, while the changes in the terrestrial water budget are displayed in Fig. 4.

a. Changes in energy budget and surface temperature
Changes in CO2SST lead to an increase in Q A at the surface [typically 3-5 W m 22 , with larger values during summertime when incoming radiation is maximum (Figs.3a,b)].This increase is caused mainly by increased incoming infrared radiation (Q Ld ) associated with the higher atmospheric CO 2 .This increased Q A is associated with a surface warming of 0.61 60.28C in North America and 0.63 60.278C in Eurasia (Figs. 3c,d), over all seasons with slightly larger values during summertime.In contrast, LULCC, which reflects either deforestation or-at a minimum-shortening of the growing season (when crops replace grasses), results in decreased Q A by a few watts per square meter over land.There is a small seasonal cycle in the reduction in Q A , particularly over Eurasia, where larger reductions are simulated due to the combined effects of snow-increased albedo and larger incoming solar radiation in summer compared to winter.The response to increased surface albedo will be  ,d).The changes in temperature are larger in spring and summer, that is, during the growing season, but the uncertainty is large in December-February (DJF), particularly in North America.
There are two remarkable features that can to be highlighted from Fig. 3: 1) In the regions considered, the CO2SST-induced warming (and increase in Q A ) and LULCC-induced cooling (and decrease in Q A ) differ much less than their global mean values.This occurs despite that changes in CO2SST lead to significant mean global annual temperature change (an average over the seven models of 0.4328C globally and of 0.6258C over land only), while the response to LULCC is negligible at that scale (an average of 20.0198C over the globe and of 20.0698C over land; see Table 5).This emphasizes the global role of CO 2 (and resulting surface ocean changes) compared to the regional significance of biogeophysical effects of LULCC.It also points to the importance of not restricting the quantification of climate change to global averages (e.g., the change in mean global annual temperature), since this hides changes of similar importance over some specific (and densely populated) regions.This issue has been addressed thoroughly by Pielke et al. (2002) and Davin et al. (2007).
2) The spread among the models is larger when the models are forced with LULCC than when they are forced with CO2SST.This results from the absence of consistent change among the various models regarding the impact of land cover type on the partitioning of Q A between Q E and Q H (section 6).When the models are forced with changes in CO2SST, it is the incoming total energy provided by the atmosphere to the LSMs that is perturbed and not necessarily the way this energy is subsequently partitioned.The increased incoming energy results in increased Q A that leads to warming in all models, since more than 75% of it is used to warm up the land.When the models' external forcing is LULCC, it is not the total incoming energy that is perturbed (except potentially via resulting feedbacks), but the way it is partitioned into absorbed and reemitted energy.Not only is this absorbed energy (and thereby Q A ) systematically reduced in all models during all seasons due to increased surface albedo but the fraction of it that is used to warm up the land is also modified by LULCC (see section 6).
In the regions of maximum LULCC changes (i.e., North America and Eurasia), the CO2SST-induced warming is compensated for by the LULCC-induced cooling, resulting in almost no change between the periods.This has strong implications for the interpretation, detection, and attribution of the observed changes between preindustrial times and present day.If the model used to detect and attribute the observed changes to various causes does not include LULCC, then erroneous conclusions may be drawn from the analyzed simulations.

b. Changes in the surface water balance
The changes in the surface water balance, defined as the difference between rainfall (P) and total evapotranspiration (E) in each model, for each season and each region, shows a distinct seasonal cycle of the simulated anomalies in response to CO2SST.The response to LULCC shows no such cycle (Fig. 4).As a result of LULCC, P 2 E increases slightly in the majority of models during winter, spring, and fall in response to a small reduction in E and no significant change in P.During summertime in North America, more than half of the models show decreased FIG. 4. As in Fig. 3, but for P 2 E (mm day 21 ).TABLE 5. Changes in mean annual global 2-m air temperature, together with the mean annual values computed over land only, in all models and the average over the seven models.''LULCC'' shows the impact of LULCC between PI times and PD, while ''CO2SST'' refers to the impact of the combined changes in atmospheric CO 2 (and equivalent GHGs) and SST and sea ice.P 2 E. However, Fig. 4 shows a very large level of spread between the models, much larger than the spread resulting from the CO2SST forcing.The reason for this increased dispersion (which was also referred to in the previous section) will be discussed in section 6.

Climate/vegetation model
When the models are forced with changes in CO2SST, a summer drying results from a combination of reduced rainfall and increased evapotranspiration.The latter is mainly a response to warmer temperatures (see, e.g., Wetherald and Manabe 2002;Gerten et al. 2007).For the summer season, for a number of models (less than half of them), the changes resulting from CO2SST are opposite of those resulting from LULCC.This again has consequences for detection and attribution efforts, since our results highlight the fundamentally different fingerprint LULCC has on the water availability in comparison to changes in CO 2 .

Changes in surface properties, fluxes, and temperature
Changes in LULCC affect many land surface characteristics (e.g., albedo, roughness length, foliage density), some of which affect Q A .Changes in these characteristics, and changes in Q A , affect the total sum of Q H and Q E , and the partitioning of these two fluxes.These changes affect surface temperature, boundary layer profiles of water and heat, and potentially cloud, convection, rainfall, etc. (Betts et al. 1996;Pitman 2003).

a. Resulting changes in land surface characteristics
In all models, in all seasons, the mean land surface albedo is higher under modern land cover than in preindustrial times (Fig. 5), especially during winter, when snow cover affects some of the temperate regions that have experienced the largest LULCC.The magnitude of the albedo change varies from less than 1% in winter in some models (e.g., ECHAM5-JSBACH) to about 8% in others (e.g., SPEEDY-LPJmL, ARPEGE-ISBA, EC-EARTH-TESSEL).The spread is quite smaller during summer (from 0.05 to about 3% at most) and fall, when vegetation foliage is fully developed and snow is not a factor.
In summer (when snow plays no role), the spread in albedo among all but one model is nearly proportional to the scale of deforestation the various models have undergone (Fig. 6).CCAM-CABLE shows a small albedo change despite a significant deforestation.In this version of the model, the parameters used in calculating canopy albedo (Sellers et al. 1992) do not vary as a function of plant functional type, making the model albedo insensitive to changes in the vegetation structure (this has been revised in more recent versions of the model).For the other models, the JJA albedo changes are roughly proportional to the deforestation scale, with an average albedo increase of 7% for an almost complete forest clearance.This proportionality is quite interesting because it means that even though the LSMs have been developed independently, the albedo's response to deforestation is quite similar from one model to another.It also implies that, to first order, modelers can perform a presimulation test of their implementation of LULCC before running experiments.One could indeed require the change in forest area to be agreed upon between modeling groups, since this largely constrains the resulting change in albedo.In DJF the interactions between vegetation and snow significantly complicate the relationship between the change in forest fraction and albedo (Fig. 6, left panel).The changes in leaf area index (LAI; Fig. 7) are not as homogeneous among the models as the changes in surface albedo.All models show decreased foliage development during wintertime in North America and Eurasia as expected, since forests and grasslands have commonly been replaced by crops, which have negligible foliage development at that time of the year.However, four out of the seven models show at least one season during which the LAI is increased once crops and pasture have replaced forests and grasses (ARPEGE-ISBA, CCAM-CABLE, ECHAM-JSBACH, IPSL-ORCHIDEE).The reasons for this increased LAI in spring (for IPSL-ORCHIDEE) or in summer (for ARPEGE-ISBA and ECHAM5-JSBACH) results either from a shift in the seasonal cycle (leaf phenology)-crop leaves in those models do not emerge at the same time as tree leaves, and the length of the growth season is quite shorter for crops than for most other plant types-or from a larger LAI for crops than for grasses (e.g., for ARPEGE-ISBA).
Regarding the relation between the scale of deforestation and the amplitude of albedo changes, there is no such clear link between the deforested fraction and changes in LAI (Fig. 8).Many variables determine the representation of LAI in most models.For example, in IPSL-ORCHIDEE, ECHAM-JSBACH, and SPEEDY-LPJmL, LAI is computed for each vegetation type and results from multiple simulated processes (photosynthesis, respiration, allocation) determined by parameters specific to individual plant function types.In all other models, LAI is prescribed per pixel and per vegetation type and is therefore independent of the potential surface climate change that the climate model may have undergone between preindustrial times and present day.

b. Changes in available energy
Following the changes in surface albedo, Q A decreases in the temperate regions, for all seasons and all models (Fig. 9).This decrease in Q A is proportional to the albedo change (not shown), and therefore approximately proportional to the scale of deforestation each model has undergone, except for CCAM-CABLE, as discussed in the previous section.In contrast to the albedo change, the changes in Q A are at a maximum during spring and summer for all models when incident solar radiation is highest.
In addition, the LULCC simulations induce a reduced atmospheric longwave radiation emission.This (indirect) effect is led by the changes in the surface energy budget, which result in a reduced heat transferred back into the upper levels and in a lower equilibrium atmospheric temperature (Van der Molen et al. 2011).
Averaged over all models, the decrease in mean annual Q A varies between 1 and more than 10 W m 22 (for EC-EARTH-TESSEL) in the northern temperate regions of North America and Eurasia, which represents a change of about 1%-10% compared to the preindustrial simulated values (Fig. 10).

c. Changes in surface heat fluxes (turbulent and thermal radiative)
The simulated decrease in Q A discussed above is accompanied-for all models, in all seasons and locationsby a systematic decrease (Figs.11a,c  Even for the five (four in Eurasia) models that decrease Q T less than Q A , the relative change (DQ T /DQ A ) varies from one model to the other.In Northern America it ranges from about 25% for ARPEGE-ISBA and ECHAM-JSBACH to more than 50% for EC-EARTH-TESSEL and CCSM-CLM and to about 100% for SPEEDY-LPJmL.Moreover, for two out of those five models (EC-EARTH-TESSEL, SPEEDY-LPJmL), the simulated change in Q T is larger than what was expected, from the simple assumption that the fraction of energy used for turbulent fluxes is the same for all periods (gray shaded area in Figs.11a,c).These results suggest that changing land cover and its associated characteristics have led to a change in the functioning of the soil-vegetation-atmosphere in a number of models.In most cases, crops and grasslands are less efficient than trees in transferring energy back to the atmosphere in the form of turbulent fluxes.This is further confirmed by Fig. 12, which shows a systematic decrease in the ratio: Q T /Q A for all seasons and for all models.Interestingly, all models at preindustrial times show rather comparable use of Q A for turbulent fluxes for all seasons, with a maximum use during summertime (;25%), followed by spring (;20%) and fall (;12%), and with minimum use during wintertime (;5%).The largest changes in the use of this energy following the imposed LULCC are obtained for EC-EARTH-TESSEL, SPEEDY-LPJmL, and IPSL-ORCHIDEE, that is, the models that undergo the largest deforestation rates (240%, 231%, 228% respectively), while the smallest changes are obtained for CCSM-CLM and ECHAM5-JSBACH, which experience the lowest deforestation rates (13% and 12%, respectively).Variations among the models can be found in the timing of the maximum change in the Q T /Q A ratio.For two models the maximum reduction of this ratio occurs in summer (IPSL-ORCHIDEE and CCAM-CABLE); for SPEEDY-LPJmL, ARPEGE-ISBA, and ECHAM5-JSBACH, the maximum changes occur in winter; while for EC-EARTH-TESSEL, it occurs in spring.
The decrease in the Q T /Q A ratio for all models at all seasons suggests that deforestation leads to an increased portion of Q A that is used to warm up the land (reduced long-wave cooling), while the sum of the turbulent fluxes (Q H and Q E ) decreases.This is a common feature shared by all models even though the season of maximum decrease varies from one model to the other.One possible cause is a decrease in the aerodynamic roughness length, which reduces the capacity of the land to exchange energy with the atmosphere via turbulence.

d. Changes in surface temperature
All models that undergo a change in their forest fraction that is larger than 15%2 simulate cooler ambient air temperature in all seasons (Fig. 13).This results from the simulated decrease in both surface albedo and Q A .One exception is IPSL-ORCHIDEE, which simulates a small warming during summer and fall in response to the large decrease in turbulent fluxes that exceeds the decrease in Q A , as discussed earlier.CCSM-CLM, ECHAM5-JSBACH, and SPEEDY-LPJmL (for Eurasia only) experience the lowest deforestation rates (between 10 and 15%) and exhibit smaller changes in temperature.Those changes though still tend to be negative.
Our results then tend to point to a rather systematic (and therefore potentially robust) cooling associated with LULCC in the temperate regions, for six out of seven models.The IPSL-ORCHIDEE model shows warming instead of cooling during summer and fall due to 1) the very strong subsequent partitioning of Q A into Q H rather than Q E in all periods (see section 5) and 2) the strong sensitivity of Q E to vegetation change through changes in surface parameters (notably, roughness length).

How do LULCC modify the partitioning of
available energy in latent and sensible heat fluxes?Is there any consistency among the models?
The partitioning of turbulent fluxes into Q E and Q H (expressed by the Bowen ratio BR 5 Q H /Q E ), and its changes since preindustrial times are shown in Fig. 14.The spread among the models is large, both in terms of its seasonal cycle and its response to deforestation.SPEEDY-LPJmL and IPSL-ORCHIDEE tend to maximize BR in all seasons except winter.From March to November, Q H is 1.1-1.5 times larger than Q E for those models, while for CCAM-CABLE, CCSM-CLM, EC-EARTH-TESSEL, and ARPEGE-ISBA, Q H is always smaller than Q E (from 0.39 to 0.71).However, this grouping of models does not imply they will similarly respond to deforestation.SPEEDY-LPJmL and IPSL-ORCHIDEE show a consistent decrease in BR in winter and increases in summer and fall, but they show opposite responses in spring.Both CCSM-CLM and EC-EARTH-TESSEL increase the return of energy in the atmosphere in terms of Q E rather than Q H , although their forest fraction has strongly diminished.For ARPEGE-ISBA, CCAM-CABLE, and ECHAM5-JSBACH, BR decreases marginally in summer, while changes in the other seasons show increases or decreases.There is therefore no consistent pattern among the various models regarding how Q E and Q H change from one season to another.Moreover, the reasons why the models differ so considerably vary from model to model.We discuss the reasons for North America, from the results displayed in Fig. 14 d The changes in CCAM-CABLE are very small.To a large degree, this is related to the small change in Q A energy (Fig. 9), which relates to the use of parameters in the calculation of albedo that do not vary as a function of vegetation type.Thus, a change in vegetation type does not cause a change in albedo and by implication Q A , and thus this element of the driver of the impact of LULCC is not captured.
d The increased Q E in CCSM-CLM arises from the partitioning of Q E into transpiration, evaporation of intercepted water, and soil evaporation.For example, summer transpiration generally decreases with the conversion of forest to crop, as does interception.Soil evaporation increases because the decrease in LAI allows more solar radiation to reach the ground.Soils in CCSM-CLM are wet, and the increased radiation produces increased evaporation.Increases in precipitation as a result of land cover change additionally increase Q E , especially in regions of North America.
d In SPEEDY-LPJmL, crops have less access to soil water, as their roots are concentrated (80%) in the upper soil layer, while trees can have more roots (up to 40%) and thus access to water in the lower soil layer.Reductions in roughness length are responsible for more stable atmospheric conditions, which also reduces Q E .During winter, the albedo increase is much larger than during summer for most models.The reduction in turbulent fluxes is more reflected on Q H than on Q E , which is always quite small in winter anyway.
d The largest signal with ARPEGE-ISBA comes from the increased snow cover over crops in winter and early spring, reducing the available energy.This reduction leads to a decrease in sensible heat flux because vegetation is mostly inactive.In summer and fall, the higher albedo of crops reduces the available energy and induces a decrease in sensible and latent heat fluxes.The reduction is larger for sensible heat flux in summer, while larger for latent heat flux in the fall.In ARPEGE-ISBA, crops have access to a smaller soil water reservoir than natural vegetation.This limitation is most obvious in the fall, when soil moisture is depleted.
d For ECHAM-JSBACH we see consistent signals in winter, spring, and summer: a decrease (increase) in LAI reduces (increase) the latent heat flux with the converse signal in the sensible heat flux.In fall the behavior is different: the almost unchanged LAI is consistent with almost no change in the latent heat flux but, surprisingly, goes along with a substantial reduction in sensible heat flux, so that maybe nonlocal effects play a role.

Discussion and conclusions
LUCID undertook a suite of climate model simulations, designed to diagnose and quantify the robust impacts of LULCC on climate between the preindustrial times and the present.Pitman et al. (2009) examined the spatial distribution of summer changes over the whole globe and concluded that there was a lack of consistency among the models.Various reasons for the discrepancies were suggested, without being able to attribute the changes to one or the other at that stage.This paper provides the detailed analysis of the LUCID results to identify the reasons for the discrepancies identified by Pitman et al. (2009).
The large spread of our model results both in terms of preindustrial and present-day values of latent and sensible heat fluxes, and their changes under deforestation confirm that the uncertainty in the Q E response is a key uncertainty in the LULCC forcing at temperate latitudes (Bonan 2008).Our results show that it is not only the change in Q E that matters but also (and probably mainly) how each individual model calculates the seasonally and surface-type-dependent Bowen ratio.
The model spread is quite surprising, given that many of the LSMs have undergone exhaustive validation of their fluxes via offline evaluation studies.However, there are a series of profound limitations in the ways LSMs are usually evaluated and intercompared.These include the following: 1) The evaluation of LSMs tends to be limited in scope, not always objective, relative to a previous version of a model (as distinct to the information content in the observed data), and rarely relative to a benchmark (Abramowitz 2005;Abramowitz et al. 2008).2) Almost all LSM evaluation is uncoupled from an atmosphere (offline); this is necessary but insufficient to establish the utility of a model, because the forcing of the atmosphere is very strong and might overshadow some differences in the behavior of the models.3) Almost all LSM evaluation has been focused on whether the mean state, or diurnal cycle, or seasonal cycle can be captured.Again, this is necessary but insufficient.Rarely is the land surface models' ability to simulate the impact, feedbacks, and associated responses to a perturbation (via elevated radiative forcing, LULCC, a severe rainfall anomaly, a drought, a heat wave, etc.) explored.However, LSMs must be able to simulate these kinds of phenomena to provide climate models the land surface influence required in climate modeling and climate projection.
Our finding therefore suggests that offline evaluation of LSMs is necessary but insufficient and should be complemented by the evaluation of the same models at least coupled to an atmospheric column.Further, an evaluation of a LSM's capacity to capture changes, as distinct from a mean state, has emerged as a priority.
The more in-depth analysis presented in this paper, compared to Pitman et al. (2009), suggests a number of robust common features shared by all models, including the following: d LULCC leads to a systematic increase in surface albedo in all seasons.For most models, this increase is proportional to the amount of deforestation imposed on the individual models.This proportionality implies that while the LSMs have been developed independently, the models respond quite similarly to LULCC in terms of albedo (7% for a full transition from forest to crop/grassland).
d This larger surface albedo causes a decrease in Q A (computed as the sum of absorbed solar energy and incident atmospheric infrared radiation); Q A decreases everywhere in the temperate regions, in all seasons and for all models.The decrease in Q A is proportional to the amount of deforestation imposed on a given model.
d The simulated decrease in Q A is accompanied, for all models, in all seasons and locations, by a systematic decrease in the sum of Q E and Q H .
d All models show a similar seasonal cycle in the amount of Q A being used for turbulent fluxes (i.e., similar Q T /Q A ).This ratio is always decreased in response to LULCC in all seasons, but this decrease varies across the model ensemble.The mean annual decrease is proportional to the intensity of LULCC.
d In most cases, crops and grasslands are less efficient than trees in transferring energy to the atmosphere in the form of turbulent fluxes due to a lower aerodynamic roughness length.
These common features and their dependence on the amount of LULCC prescribed in each model suggest that, for an agreed amount of deforestation that occurred over specific periods, the dispersion among the models will be significantly smaller.However, some persistent disparities remain.First, there is no consistency in how Q T is partitioned between Q E and Q H throughout the annual cycle.Second, there is no consistency in the seasonal response of Q H and Q E to LULCC.In some LSMs, deforestation leads to a decrease in both fluxes in response to the decrease in Q A and Q T , while in others Q H increases and Q E decreases (or vice versa).This is dependent on how these processes are represented in each LSM.We are not able to attribute these differences to the many parameterization variations.However, our conclusions do point to an urgent need to revisit the way LSMs are evaluated if we are to resolve why they disagree on the impact of LULCC.One significant problem we are facing is that LSMs are generally evaluated offline, forced with prescribed atmospheric forcing.Our results suggest that this is necessary but insufficient.Evaluation of LSMs has to account for atmospheric feedbacks (e.g., Santanello et al. 2009).An evaluation framework should then be completed by series of analyses that determine how well LSMs capture the contrasting dynamic properties of various vegetation types, which are relevant for biosphere-atmosphere interactions (e.g., water use efficiency, dynamics of evaporative fraction, effective temperature sensitivity of carbon balance).
Increased concentration of greenhouse gases in the atmosphere, and the subsequent changes in sea surface temperatures and sea ice extent, are often used as the main drivers of climate change also over land.Our results suggest that such an assumption leads to erroneous conclusions regarding the land surface impacts of climate change in regions where LULCC has been significant.LULCC affects a number of variables to a similar magnitude, but of opposite sign, in increasing greenhouse gas concentrations.LULCC therefore has the potential to mask a regional warming signal, with the resulting risk that detection and attribution studies may miss a clear greenhouse signal or misattribute a greenhouse signal if LULCC is poorly accounted for.Detection and attribution is a complex process that is beyond the scope of this paper (see, e.g., Stott et al. 2010).However, our results suggest that including LULCC could improve the regional-scale detection of the impacts of specific forcings by ensuring that land cover's contribution to any regional changes is appropriately represented.LULCC will suppress the impacts of, for example, increasing CO 2 in some regions that cool due to land cover change and amplify the impacts of increasing CO 2 in regions that warm due to land cover change.In the former case, there is a risk of missing the detection of a CO 2 signal, while in the latter there is a risk of a false-positive detection of a CO 2 signal.Aerosols, which typically cool, particularly strongly at regional scales, are an additional forcing that might be masked by the misrepresentation of LULCC.
Our findings argue for the inclusion of LULCC in climate projections, as now in process for the CMIP5 simulations (e.g., Arora et al. 2011).However, we have also shown that the differences among the seven models' response to LULCC is larger than the differences that results from the change in CO2SST.Since LULCC is implemented in most CMIP5 models, we expect a larger divergence among climate models in comparison to earlier efforts over regions of intense LULCC, coincident with dense human populations.This problem will only be resolved via a more systematic effort within the climate modeling and land surface modeling communities, starting with a coordinated evaluation of how to represent LULCC, how well models capture the impacts of LULCC in both offline and coupled simulations.

FIG. 1 .
FIG.1.Changes in the extent covered with crops and pasture between PD and PI times.Yellow and red are used when the extent of anthropogenic areas have increased since PI times, while blue refers to abandoned lands.The two contours that are drawn on the map highlight the regions that will further be used for specific analysis (North America and Eurasia).

FIG. 2 .
FIG. 2. (a)The 1870 extent (in fraction of total area) covered by crops (gray), grassland types (orange), evergreen trees (green), deciduous trees (blue), and desert (white) for all seven models and two different geographical locations (North America and Eurasia; located in Fig.1).(b) Differences (in fraction of total area) in each of those vegetation types between PD and PI times.The dashed black line on both graphs shows the fractional coverage by crops that have been provided.

FIG. 3 .
FIG. 3. Box-and-whisker plots of the simulated changes between the PI period and PD in (a),(b) available energy (W m 22 ) and (c),(d) surface air temperature (8C) for all seasons and for (a),(c) North America and (b),(d) Eurasia.The mean ensemble values of each individual model and each set of experiment (i.e., PD 2 PIv and PDv 2 PI for the CO2SST impacts; PD 2 PDv and PIv 2 PI for the LULCC impacts) have been used to create this plot.The bottom and top of the box are the 25th and 75th percentiles, and the horizontal line within each box is the 50th percentile (the median).The whiskers (straight lines) indicate the ensemble maximum and minimum values.

FIG. 5 .
FIG.5.Seasonal LULCC-induced changes in the simulated surface albedo (%) between PD and PI times by the models listed at the top of each panel.The shading refer to the differences calculated between the simulations that are forced with PD SSTs, CO 2 , and GHG (black, simulation PD minus simulation PDv) and PI SSTs, CO 2 , and GHG (gray, simulation PIv minus simulation PI).Presenting both dark and gray bars demonstrates the robustness of the impacts of LULCC changes, largely independent of the state of the background GHG and surface ocean's temperatures.(top) North America and (bottom) Eurasia.All seasons are plotted from (left) winter (DJF) to (right) fall [September-November (SON)].

FIG. 6 .
FIG.6.LULCC-induced changes in the simulated surface albedo (%) plotted against the changes in forest fraction between PD and PI times.Symbols refer to the various models used (ARP: ARPEGE-ISBA; CCA: CCAM-CABLE; CCS: CSSM-CLM; ECE: EC-EARTH-TESSEL; ECH: ECHAM5-JSBACH; IPS: IPSL-ORCHIDEE; SPE: SPEEDY-LPJmL).Red symbols are for Eurasia, while black symbols are for North America.(left) Winter changes (DJF); (right) summer changes (JJA).Albedo values have been computed, for each model, as the mean over not only the various ensemble members but also the two sets of simulations (i.e., PD 2 PDv and PIv 2 PI).
FIG. 8.As in Fig.6, but for the simulated or imposed changes in LAI.

FIG. 10 .
FIG. 10.Mean annual LULCC-induced change in available energy, averaged over all models and all experiments [0.5 (PD 2 PDv 1 PIv 2 PI)].Shadings refer to the changes in W m 22 , while isolines represent the standard deviation between the models.

FIG. 11 .
FIG. 11.Mean summer LULCC-induced changes (W m 22 ) in (a),(c) Q T plotted against the mean summer changes in Q A ; (b),(d) longwave radiation emitted by the surface (Q LU ) plotted against the changes in the difference between Q A and Q T .(a),(b) North America; (c),(d) Eurasia.Symbols refer to individual models.Plain line in all panels represents the 1:1 regression.Gray shaded area in (left) refers to the changes in Q T that would have occurred if Q T /Q A ratios were the same at PD and PI periods (see section 4c for discussion).

FIG. 12 .
FIG. 12. Mean seasonal values of (a) the fraction of available energy used for turbulent fluxes (i.e., Q T /Q A ; %) for each LSM for PI conditions.(b) The LULCC-induced changes in this fraction between PD and PI periods.

FIG. 14 .
FIG. 14. Mean seasonal values of (a) the BR for each LSM for PI conditions.(b) The LULCC-induced changes in this fraction between PD and PI periods.

TABLE 1 .
List of climate (or atmosphere only) models and associated LSMs used in the first LUCID set of experiments.
(CABLE)to more than 30% in the earth system model based on European Centre for Medium-Range Weather Forecasts (ECMWF) modeling systems (EC-EARTH)-Tiled ECMWF Scheme of Surface Exchanges over Land

TABLE 3 .
Strategy used in the various models to implement the crops and pasture changes onto their background land cover map.CRU TS 2.0 5

TABLE 3
or pastures) in the remainder of the grid cell (e.g., SPEEDY-LPJmL).

TABLE 4 .
Forest extent (km 2 ; numbers in parenthesis represent the fractional area covered by forests) at both PI time and PD, together with changes in the forested areas (km 2 and %) between those two time slices.All models are shown for both geographical regions.The magnitude of the simulated changes is generally smaller over Eurasia than over North America due to smaller LULCC changes in the former region (see Table4).The dominant impact of lowered Q A is a quasi-systematic cooling of the land by about 20.448 60.48C in North America and about 20.3 60.38C in Eurasia in all seasons(Figs.3c The seasonal cycle of leaf area index in IPSL-ORCHIDEE is slightly shifted toward earlier dates (section 5a), which leads to larger values of LAI in spring and smaller values in all other seasons.ORCHIDEE is sensitive to changes in LAI, and Q E is increased in spring (despite the decrease in Q A ) and decreased in all other seasons.The warming of the ambient air in summer and fall results from this significant reduction in Q E that more than compensates the change in Q H .In June-August (JJA), EC-EARTH mimics the effects of ORCHIDEE in spring: the reduction in Q A gives rise to excessive reductions in Q H and a compensating increase in Q E , supported by lower evaporation and soil moisture depletion rates in DJF and March- : d d